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Above tens of GeV, 7-ray observations with the Fermi Large Area Telescope (LAT) can be dominated by 
statistical uncertainties due to the low flux of sources and the limited acceptance. We are developing a new 
event class which can improve the acceptance: the “Calorimeter-only (CalOnly)” event class. The LAT has three 
detectors: the tracker, the calorimeter, and the anti-coincidence detector. While the conventional event classes 
require information from the tracker, the CalOnly event class is meant to be used when there is no usable 
tracker information. Although CalOnly events have poor angular resolution and a worse signal/background 
separation compared to those LAT events with usable tracker information, they can increase the instrument 
acceptance above few tens of GeV, where the performance of Fermi-LAT is limited by low photon statistics. In 
these proceedings we explain the concept and report some preliminary characteristics of this novel analysis. 


1. Introduction 

The Fermi Large Area Telescope {Fermi LAT ) 
is an instrument on the Fermi 7 -ray telescope oper¬ 
ating from 20 MeV to over 300 GeV. The instrument 
is a 4 X 4 array of identical towers, each one consist¬ 
ing of a tracker-converter (TKR), based on Silicon 
detector layers interleaved with Tungsten foils, where 
the photons have a high probability of converting to 
pairs, which are tracked to allow reconstruction of the 
7 -ray direction and a segmented calorimeter (CAL), 
made of Csl crystal bars, where the electromagnetic 
shower is partially absorbed to measure the 7 -ray en¬ 
ergy. The tracker is covered with an anti-coincidence 
detector (ACD) to reject the charged-particle back¬ 
ground. Further details on the LAT, its performance, 
and calibration are given by [T] and [ 2 ]. 

Most of the science done with Fermi LAT spans 
photons with energies from 50 MeV to about 10 GeV, 
where the sensitivity of the instrument is good and the 
available number of detected photons high. However, 
there are many sources which emit 7 -rays above a few 
tens of GeV. These energies that are almost accessi¬ 
ble by the current generation of Imaging Atmospheric 
Cherenkov Telescopes (lACTs). Even though the de¬ 
tection area of LAT is small (in comparison to that of 
lACTs), Fermi LAT provides all-sky coverage and a 
very high duty cycle, which are crucial characteristics 
for producing 7 -ray source catalogs and study source 
variability in an unbiased way. A prime example is 
the the first Fermi LAT catalog of >10 GeV sources 
(IFHL) m, which contains 514 sources, out of which 
~100 sources have already been detected at very high 
energies(>100 GeV, or VHE), and ^200 additional 
sources have been identified as good candidates to be 
VHE emitters and be detected with lACTs. 
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Figure 1: Schematic diagram of the Large Area 
Telescope. The telescopes dimensions are 
1.8m X 1.8m x 0.72m. [T] 


The performance of Fermi LAT above 10 GeV 
is excellent. The angular resolution and sig¬ 
nal/background separation is best at the highest pho¬ 
ton energies, where one only suffers from a slight dete¬ 
rioration of the energy resolution due to the fact that 
the showers are no longer contained in the calorimeter. 
However, the steep falling photon flux with energy of 
most 7 -ray sources, together with the relatively small 
effective area of LAT (^Im^), results in a substan¬ 
tial limitation due to the very low number of detected 
photons (e.g., in the IFHL, many sources were char¬ 
acterized with only 4-5 photon events over a back¬ 
ground signal of 0-1 event). The low statistics from 
7 -ray sources is going to be an even larger problem for 
the second Fermi high-energy LAT catalog (2FHL, 
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in preparation), which is expected to consist on 7 -ray 
sources detected above 50 GeV (instead of 10 GeV). 

In these proceedings we report an analysis which 
can help increase the photon statistics at few tens of 
GeV, hence improving the ability to perform science 
at the highest LAT energies, where the lAGTs start 
operating. The methodology is still being developed. 
Here we only present the concept and report some 
preliminary characteristics. 


2. The Calorimeter-only (CalOnly) 
Fermi-LAT analysis 

The regular Fermi LAT event classes require us¬ 
able information from the TKR. This is a sensible 
approach, given that the TKR information is crucial 
to determine accurately the incoming direction of the 
7 -ray event. The LAT TKR comprises only ^1.5 ra¬ 
diation lengths (on axis), which means that a large 
fraction of 7 -rays from the astrophysical sources are 
discarded at the very beginning of the analysis because 
they do not convert in the TKR, or they convert in 
the bottom layers and the TKR information is not 
sufficient for a proper determination of the incoming 
direction of the 7 -ray event. This situation is depicted 
in Figure 

In the standard LAT analysis, the CAL is essen¬ 
tially used for signal/background separation (together 
with the TKR and ACD) and to determine the en¬ 
ergy of the 7 -rays and electrons. The LAT has a ho- 
doscopic calorimeter, consisting of 16 towers with 8 
layers of 12 crystals each because of three driving rea¬ 
sons. First, shower profiling improves energy recon¬ 
struction. Next, shower topology contains valuable 
information for signal/background separation. Last, 
it realizes independent event acceptance and recon¬ 
struction. 

To create a usable event class without using TKR 
information, one must determine the incoming direc¬ 
tion of the 7 -rays with sufficient resolution (a few de¬ 
grees) while keeping a reasonable background rejec¬ 
tion (^0.999).This would increase the number of avail¬ 
able high-energy events for performing 7 -ray astron¬ 
omy. The above mentioned (non-standard) LAT anal¬ 
ysis, dubbed Calorimeter-only (CalOnly) analysis, is 
currently being developed by the Fermi LAT collab¬ 
oration, and is aimed to deliver yet another class of 
events, the CalOnly event class, which may be added 
to the other photon event classes coming from the reg¬ 
ular LAT analysis. For the CalOnly analysis, it fully 
reconstructs the electromagnetic showers and deter¬ 
mine its main axis, which points to the direction of 
the incoming 7 -ray event. 

Since the main event trigger for LAT is based on 
the TKR, those events with no usable TKR informa¬ 
tion have a low chance of being recorded and trans- 
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Figure 2: The left panel depicts two events that convert 
to electron pairs in the TKR (from the broken lines to 
the solid lines in these figures). These events would be 
used in the regular Fermi LAT analysis. The right panel 
shows two events without usable TKR information. One 
is a side-entering event that goes through a small fraction 
of the TKR. The other one crosses the entire TKR, but 
converts to an electron pair in the CAL. These two 
events would be rejected and not be used for regular 
LAT analysis. These are the type of events which could 
be recovered with a dedicated analysis that does not 
require TKR information, the so-called CalOnly events. 


ferred to the ground. However, large energy deposi¬ 
tions in CAL generate a trigger that is fully indepen¬ 
dent of TKR, and on-board event selection records 
events that deposit energy larger than 20 GeV in CAL. 
Consequently, the CalOnly event class will only be ef¬ 
fective above a few tens of GeV. 


2.1. Pass 8 

The development of the CalOnly analysis is done in 
the context of Pass 8 , which is the new iteration of the 
LAT event-level analysis package. Pass 8 was orig¬ 
inally designed to address the effect of coincidences 
with cosmic rays (‘ghost’ events), but quickly evolved 
in a comprehensive revision of the instrument simu¬ 
lation, the event reconstruction, and the background 
rejection, with the goal of improving all the aspects of 
the LAT performance: larger acceptance, better angu¬ 
lar and energy resolution, and extension of the energy 
reach below 100 MeV and in the TeV range. 

The details on the Pass 8 analysis chain can be 
found in [4]. The CAL reconstruction begins with a 
clustering stage that tries to isolate the genuine 7 -ray 
shower from smaller energy deposition due to ghost 
events. At this point we can exploit the segmentation 
of the CAL to identify the energy deposition centroid 
and the shower axis (via a moment analysis) that, 
for CalOnly, corresponds to the photon incoming di¬ 
rection. This direction is propagated to the ACD (in 
addition to the tracks from the tracker) in order to as¬ 
sociate energy deposition in the tiles and discriminate 
charged particle without the TKR direction. This is 
one of the main improvements introduced with Pass 
8 and proved to be very useful even if its separation 
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log10(Combined energy and position deviation in ACD) 


Figure 3: The left, middle and right panels show the logarithm of the CAL first cluster transverse profile RMS, the 
logarithm of of fhe profile fit computed over a 100 mm radius cylinder around the trajectory, and the logarithm of 
combined energy and position deviation in ACD respectably. The last one needs some more description. It is the 
number of sigmas less than an expected MIP signal, combined with the number of sigmas the track propagation is 
away from tile or ribbon most likely to veto the first CAL cluster and the number of sigmas the energy deposited in the 
ACDs is away from the expected amount. If the value of this combined quantity is zero, then the log 10 of this quantity 
is set to -4. This is more likely to happen to MIPs than for 7 -rays. If there are no tracks associated to ACD signals, 
this quantity is set to +4. Charge particles are likely to have smaller values than the 7 -rays. 


capability is limited by the angular resolution of the 
CAL, that is obviously worse than that of the TKR.. 

Another important improvement in the CAL re¬ 
construction is the energy measurement that, for 
CalOnly, is based on a full three-dimensional fit of 
the shower energy deposition. This method needs a 
precise modeling of the longitudinal and lateral de¬ 
velopment of showers inside the CAL and a reference 
axis. The latter, usually taken from the tracker, can 
be obtained from the aforementioned moment analysis 
with a small change in performance. 

The last step of the Pass 8 development is the high- 
level analysis that links together all the outputs of 
the reconstruction and classify events as good 7 rays 
or not. This is the core of the CalOnly development 
and is described in next section. 

2.2. Signal / Background separation 

The Fermi LAT needs to reject a cosmic-ray back¬ 
ground that outnumbers the signal (7 rays) by many 
orders of magnitude, and hence an efficient sig¬ 
nal/background separation (rejecting 10 ^“^ of the 
background events) is required to be able to per¬ 
form 7 -ray astronomy. The LAT background consists 
mostly on protons and electrons, but also on alpha 
particles and heavy nuclei. As the 7 -ray energy in¬ 
creases, we have a natural improvement in the sig¬ 
nal/background ratio due to the fact that most 7 -ray 
sources have spectra that can be parameterized with 
power-law indices harder than 2.5 (often even harder 
than 2 . 0 ), while the spectra of the proton background 
follows a power-law index of ^2.7 and that of the elec¬ 
trons a power-law index of ^3.1. 


The rejection of the background cosmic-ray events 
in the CalOnly analysis is based on the different topol¬ 
ogy of electromagnetic and hadronic showers, and the 
ACD signals produced by the charged particles. It 
is worth noting that, while protons and heavy nuclei 
can be effectively distinguished from 7 rays using only 
information from the CAL, the electrons/positrons 
produce electromagnetic showers that are essentially 
identical to those of the 7 rays, and hence the informa¬ 
tion from the ACD is crucial to be able to reject elec¬ 
trons. The left and middle panels in Figures show 
the normalized MC distributions of two CAL-related 
parameters that can effectively distinguish between 
electromagnetic {gammas and leptons) and hadronic 
showers. In order to be able to reject the leptons, one 
needs the help of ACD-related parameters, as depicted 
in the right panel in Figure 

In order to maximize the separation of signal and 
background, instead of making simple cuts in the dis¬ 
tributions of CAL and ACD parameters as the ones 
shown above, we perform the analysis through a multi¬ 
variate analysis (MVA) that uses a large number of 
CAL and ACD parameters. For the most effective 
background rejection, we applied the Boosted Deci¬ 
sion Tree (BDT), one of the methods of multi-variate 
analysis. In this method, we train many classifica¬ 
tion trees with Monte Carlo (MC) data, for which 
can identify unambiguously what is signal and what 
is background. We are using the ROOT-based TMVA 
package to train the classification tree analysis [5] . We 
build the trees and then evaluate the gamma-likeness 
of each event. Next we can cut on the gamma-likeness 
and get the events classified as signal or background. 
By selecting events with a very high gamma-like ness. 
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Figure 4: Efficiency of MC 7 events (blue) and rejection 
of MC background events (red) vs. cut value of a 
classifier output. 

one can increase the purity of the selected data set, 
but at the expense of reducing the number of 7 -ray 
candidates. This is represented in FigureThe op¬ 
timal value to increase the signal/background is typ¬ 
ically obtained for a cut value between 2 and 3 (de¬ 
pendent on the energy range and incidence angle of 
the 7 -ray considered). 

At the present time we are optimizing the classi¬ 
fication tree by modifying the input CAL and ACD 
parameters (including creating new composite vari¬ 
ables), as well as by building the classification trees in 
different modes. Consequently, the results presented 
in this proceedings should be considered as prelimi¬ 
nary, and likely reporting a lower limit of the actual 
performance of the CalOnly analysis. 

2.3. Quality of the reconstructed events 

In this section we address the quality of the re¬ 
constructed CalOnly events (after signal/background 
separation) using dedicated MC simulations of 7 -ray 
events. 

Two basic quantities are being evaluated: the an¬ 
gular and the energy resolution. Given that the thick¬ 
ness of the calorimeter increases rapidly with the in¬ 
cidence angle, one expects a different performance for 
low and high incidence angle 7 rays. In this section 
we define low (high) incidence angle as smaller (larger) 
than 53 degrees (cos(53°) ^ 0.6) and evaluate the per¬ 
formance for these two cases. And naturally, as it 
occurs in the regular Fermi LAT analysis, the perfor¬ 
mance can also vary with the energy of the incoming 
7 ray. Here we define low (high) energy as being in 
the range ^30-100 GeV (^100-300 GeV), and evalu¬ 
ate the performance for these two energy bands. 

Figure shows the normalized distributions in the 
error of the reconstructed directions for low/high en¬ 
ergy bands and incidence angles. The angular resolu¬ 
tion can be defined as the 68 % containment in those 


distributions (PSF 68 ), which would result in ^2 de¬ 
grees for high inclination 7 rays, and ^3-4 degrees for 
low inclination 7 rays (with a relatively small depen¬ 
dence on the energy). The PSF 68 for regular LAT 
photons (i.e. with usable TKR information) at these 
energies is ^ 0 . 1 - 0.2 degrees, which is more than one 
order of magnitude better than for GalOnly photons. 


Angular resolution 



Figure 5: Normalized histograms of the angular distances 
between reconstructed and MC direction of CalOnly 
events in 2 inclination angle x 2 energy bins. 


Energy resolution 



Figure 6 : Normalized histograms of the dispersions in 
the reconstructed energy of MC CalOnly events in 2 
inclination angles x 2 energy bins. 

Figure 1^ shows the normalized distributions in the 
error of the reconstructed energies for low/high energy 
bands and incidence angles. The energy resolution 
can also be defined as the 68 % containment on these 
distributions, and using the largest distance from the 
peak position to the edge of the 68 % containment. 
The energy resolution is ^ 3 — 4% for high inclination 
and ^ 10 — 15% for low inclination 7 -rays, with little 
dependence on the energy of the event. This perfor¬ 
mance is very close to that of regular LAT photons. 
The quality of the energy reconstruction is mainly 
connected to the path length of shower axis (related 
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to the shower containment) and the accuracy of the 
shower direction reconstruction. While the latter is 
worse for CalOnly events, this class can benefit from 
a larger field of view and therefore longer trajecto¬ 
ries. It must be noted that both direction and energy 
resolution can be improved with a dedicated selection 
of good quality events, at the price of a lower effec¬ 
tive area. The best trade off between these conflicting 
requirements is still to be evaluated. 

3. Conclusions 

Pass 8 provides an unprecedented framework to de¬ 
velop an analysis that uses events without usable TKR 
information. 

The CalOnly event class, currently under devel¬ 
opment within the Fermi-LAT collaboration, could 
be used to increase the acceptance of Fermi LAT 
above few tens of GeV (where the performance is pho¬ 
ton statistics limited), by recovering for astronomical 
studies 7 -ray events without usable TKR information. 
This implies that CalOnly events will have a worse 
signal/background separation and angular resolution, 
when compared to the regular LAT events. However, 
they can have a better energy resolution, if consider¬ 
ing the high incidence angle events. 

The CalOnly event class may be particularly rele¬ 
vant in the following two scientific topics: 

• Search for line-signals potentially coming from 
Dark Matter annihilation (because of the larger 
number of events and the excellent energy reso¬ 
lution for the large-incident angle events) 

• Study of transient events like GRBs and AGN 
flares (because of the larger number of events 
and the valuable increase in the temporal cover¬ 
age of the source) 
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